Why Ni 3 Al is an itinerant ferromagnet but Ni 3 Ga is not 
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N13AI and NisGa are closely related materials on opposite sides of a ferromagnetic quantum 
critical point. The Stoner factor of Ni is virtually the same in both compounds and the density of 
states is larger in NisGa. So, according to the Stoner theory, it should be more magnetic, and, in 
LDA calculations, it is. However, experimentally, it is a paramagnet, while M3AI is an itinerant 
ferromagnet. We show that the critical spin fluctuations are stronger than in NiaGa, due to a 
weaker q-dependence of the susceptibility, and this effect is strong enough to reverse the trend. 
The approach combines LDA calculations with the Landau theory and the fluctuation-dissipation 
theorem using the same momentum cut-off for both compounds. The calculations provide evidence 
for strong, beyond LDA, spin fluctuations associated with the critical point in both materials, but 
stronger in NiaGa than in NisAl. 
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Recent low temperature experiments on clean mate- 
rials near ferromagnetic quantum critical points (QCP) 
have revealed a remarkable range of unusual properties, 
including non-Fermi liquid scalings over a large phase 
space, unusual transport, and novel quantum ground 
states, particularly coexisting ferromagnetism and su- 
perconductivity in some materials. Although critical- 
ity usually implies a certain universality, present ex- 
periments show considerable material dependent aspects 
that are not well understood, [1] e.g. the differences 
between UGe 2 and URhGe [2,3] and ZrZn 2 , [4] which 
both show coexisting ferromagnetism and superconduc- 
tivity but very different phase diagrams, in contrast to 
MnSi, where very clean samples show no hint of super- 
conductivity around the QCP. [5] Generally, approaches 
based on density functional theory (DFT) are success- 
ful in accounting for material dependence in cases where 
sufficiently accurate approximations exist. To proceed 
in this direction it is useful to study benchmark systems 
for which detailed experimental data are available and 
which pose challenges to theory. Here we report a study 
of the closely related compounds Ni 3 Al and Ni 3 Ga. Both 
of these have the ideal cubic CU3AU cP4 structure, with 
very similar lattice constants, a = 3.568 A and a = 3.576 
A, respectively, and have been extensively studied by var- 
ious experimental techniques. N13AI is a weak itinerant 
ferromagnet, T c — 41.5 K and magnetization, M=0.23 

/cell (0.077 /Ltg/Ni atom) [6] with a QCP under pres- 
sure at P c =8.1 GPa, [7] while NisGa is a strongly renor- 
malized paramagnet [8]. Further, it was recently reported 
that M3AI shows non- Fermi liquid transport over a large 
range of P and T range down to very low T. [9] 

DFT is an exact ground state theory, and as such 
should properly describe the magnetic ground states of 
metals. However, common approximations to DFT, such 
as the LDA and GGA, are based on the properties of 
the uniform electron gas at densities that occur in solids. 



At these densities it is rather stiff with respect to spin 
fluctuations and is not close to any magnetic instabil- 
ity. As a result, the LDA description of magnetism is 
at a quasi-classical mean field level, {i.e. Stoner level), 
and neglects fluctuations due to soft magnetic degrees of 
freedom. This leads to misplacement of QCPs and over- 
estimates of the magnetic tendencies of materials near 
QCPs, as well as such known problems as the incorrect 
description of singlet states in molecules with magnetic 
ions. In fact, practically all cases where the LDA sub- 
stantially overestimates the tendency towards magnetism 
are materials near a QCP [10-13] - a fact that can po- 
tentially be used as a screen for materials with large fluc- 
tuation effects. [14] Previous LDA calculations showed 
that the magnetic tendency of both materials is overes- 
timated within the LDA, and that NiaGa is incorrectly 
predicted to be a ferromagnet. [15-20] Moreover, as our 
present results show, in the LDA the tendency to mag- 
netism is stronger in NiaGa than M3AI, opposite to the 
experimental trend. This poses an additional challenge 
to any theory striving to describe the material depen- 
dent aspects of quantum criticality. The two materials 
are expected to be very similar electronically (our results 
confirm this, and identify the small difference between 
the two as due to relativistic effects associated with Ga 
in NiaGa). Thus they offer a very useful and sensitive 
benchmark for theoretical approaches. We use this to 
test an approach based on the fluctuation dissipation the- 
orem applied to the LDA band structures with an ansatz 
for the cut-off q c . We find that this approach corrects 
the ordering of the magnetic tendencies of the materials, 
and gives the right ground states at ambient pressure as 
well as a reasonable value of P c for Ni 3 Al. 

Our LDA calculations were done using the gen- 
eral potential linearized augmented-plane-wave (LAPW) 
method with local orbital extensions [21,22] in two imple- 
mentations [22-24], with the exchange-correlation func- 
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tional of Hcdin-Lundqvist with the von Barth-Hedin spin 
scaling [25,26]. Up to 816 inequivalent k-points were used 
in the self-consistent calculations, with an LAPW basis 
set defined by the cut-off RsK max =9, plus local orbitals 
to relax linearization errors. Larger numbers of k-points 
between 2300 and 4060 were used in the Fermi surface in- 
tegrations. The LDA electronic structure is given in Fig. 
1 and Table 1, while results of fixed spin moment cal- 
culations of the magnetic properties at the experimental 
lattice parameters and under hydrostatic compression are 
given in Figs. 2 and 3. The two compounds are very sim- 
ilar in both electronic and magnetic properties, the main 
apparent difference being the higher equilibrium moment 
of Ni 3 Ga (0.79 ^s/f.u. vs. 0.71 /is/f-u.), in agreement 
with other full potential calculations. [19,20] 
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FIG. 1. Calculated LDA band structure (top) and density 
of states (bottom) per f.u. for non-spin-polarized N13AI (solid 
lines) and NisGa (dotted lines). Ef is at eV. 

The propensity towards magnetism may be described 
in terms of the Stoncr criterion, IN(Ep), where I 
is the so-called Stoner parameter, which derives from 
Hund's rule coupling on the atoms. For finite mag- 
netizations, the so-called extended Stoncr model [27], 
states that to the second order in the spin density the 
magnetic stabilization energy is expressed as AE = 
M 2 [/ M m dm/2N(m) - 7/4], where N(M) is the den- 
sity of states averaged over the exchange splitting corre- 



sponding to the magnetization M. Fitting our fixed spin 
moment results to this expression, we find Iai = 0.385 
eV and I Ga = 0.363 eV. These gives IN(E F ) =1.21 
and IN(E F ) = 1.25 for Ni 3 Al and Ni 3 Ga, respectively. 
Both numbers are larger than 1, corresponding to a fer- 
romagnetic instability, and the value for Ni 3 Ga is larger 
than that for Ni 3 Al. Importantly, the difference comes 
from the density of states, since Iai > lea- In both com- 
pounds, magnetism is suppressed by compression, with 
an LDA critical point at a value 5a/ a ~ -0.05 - -0.06. In 
Ni 3 Al, the critical point at 5a/ a =-0.058 corresponds to 
the pressure of P c =50 GPa, [28] which is much higher 
than the experimental value. It is interesting that, as 
in ZrZn2 [11], the exchange splitting is very strongly k- 
dependent; for instance, in Ni 3 Al at some points it is as 
small as 40 meV//xs near the Fermi level, while at the 
others (of pure Ni d character) it is close to 220 meV/^s. 
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FIG. 2. Energy vs. fixed spin moment for N13AI and Ni^Ga 
at the experimental lattice parameters. The energy zero is set 
to the non-spin-polarized value. 

Notwithstanding the general similarity of the two com- 
pounds, there is one important difference near the Fermi 
level, specifically, the light band crossing the Fermi level 
in the middle of the T-M or T-X directions is steeper in 
Ni 3 Al (Fig. 1). This, in turn, leads to smaller density of 
states. This comes from a different position of the top of 
this band at the T point, 0.56 eV in Ni 3 Ga and 0.85 eV 
in N13AI. The corresponding electronic state is a mixture 
of Ni p and Al (Ga) p states, and is the only state near 
the Fermi level with substantial Al (Ga) content. Due 
to relativistic effects, the Ga p level is lower than the Al 
p level and this leads to the difference in the position of 
the corresponding hybridized state. Note that this is a 
purely scalar relativistic effect. We checked that includ- 
ing spin orbit docs not produce any further discernible 
difference. 

Returning to magnetism, the fixed spin moment calcu- 
lations provide the energy E as a function of the magne- 
tization M (Fig. 2). One can write a Landau expansion 
for E{M) as 
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E(M) = a 2 M 2 /2 + ai M 4 /4: + a 6 M 6 /6 + ■■■ (1) 

Treating this as a mean-field expression and adding the 
effects of spin fluctuations [29] leads to renormalization of 
the expansion coefficients. The renormalized coefficients 
a,i are written as power series in the averaged square of 
the magnetic moment fluctuations beyond the LDA, £ 2 , 

n+i — 1 

~a2n = Y. C n+l-l a 2(n +i) e J[ (1 + 2fc/3). (2) 
i>0 k=n 

£ may be estimated by requiring that the corrected 
Landau functional (2) reproduces the experimental mag- 
netic moment (for Ni 3 Al) or experimental magnetic sus- 
ceptibility (for Ni 3 Ga). The "experimental" £'s obtained 
in this manner are are 0.47 and 0.55, respectively, which 
implies that spin fluctuation effects must be stronger in 
Ni 3 Ga than in Ni 3 Al. 

We shall now make a link between this fact and the 
electronic structures. A standard formula for estimating 
£ 2 comes from the fluctuation-dissipation theorem [30], 
which establishes that 

= (2ft/fi) J d?q I (dw/27r)Im X (q,a;), (3) 

where f2 is the Brillouin zone volume, and x is the mag- 
netic susceptibility. Using the lowest order expansion for 

X, 

Xo(q, u) = N(E F ) - aq 2 + ibuo/q (4) 

x _1 (q,w) = Xo 1 ^^) - J > ( 5 ) 

where xo (q 5 is the non-interacting susceptibility, one 
can derive a formula for £ 2 [29,30], whose coefficients can 
be related to the characteristics of the electronic struc- 
ture [30,31]. The final results reads 

e = bv ^f n F)2 [Q i Hi + Q- 4 ) + Hi + Q% (6) 

where a = (d 2 {N{E F )v 2 x ) / 'dE%) / '12, b 
(N(E F )v~ 1 )/2, v F = ^Z{d 2 {N{E F )v 2 x ), Q = q c ^/ajb^, 
and q c is the cutoff parameter for integration in Eq.3. 
The physical meaning of these parameters is as follows, a 
defines the rate at which the static susceptibility x{Q-> 0) 
falls away from the zone center, i.e. the extent to which 
the tendency to ferromagnetism is stronger than that 
to antifcrromagnetism. This translates into the phase 
space in the Brillouin zone where the spin fluctuations 
arc important, b controls the dynamic effects in spin 
susceptibility. The cutoff parameter q c is the least well 
defined quantity in this formalism. One obvious choice 
is q c = y / N(E F )/a, because for larger q the approxima- 
tion (4) gives unphysical negative values for the static 
susceptibility. On the other hand, this choice leads to 
noticeably different cutoffs for the two compounds, while 



one may argue that q c should reflect mainly the geometry 
of the Fermi surface and thus be practically the same in 
these two cases. Furthermore, the fermiology of these 
compounds is very complicated: in the paramagnetic 
state, there are 4 Fermi surfaces, two small and two large 
(one open and one closed). In this situation, it is hardly 
possible to justify any simple prescription for q c . There- 
fore, we have chosen a different route: we assume that 
q c is the same for both materials, and choose a number 
which yields a good description of both the equilibrium 
moment in N13AI and the paramagnetic susceptibility in 
Ni 3 Ga, q c = 0.382 1 . Note that this is larger that the 
diameters of the small Fermi surfaces but smaller than 
the radius of the Brillouin zone, ~ 0.5 . 

To calculate the above quantities, especially a, we need 
accurate values of the velocities on a fine mesh. Nu- 
merical differentiation of energies within the tetrahedron 
method proved to be too noisy. Therefore we use the 
velocities obtained analytically as matrix elements of the 
momentum operator, computed within the optic program 
of the WIEN package. A bootstrap method [32], as de- 
scribed in Ref. [31], was used to obtain stable values 
for a, b. We found for N13AI, using as the energy unit 
Ry, the length unit Bohr, and the velocity unit Ry-Bohr, 
a = 230, b = 210, v F = 0.20, and £ = 0.445 fi B - For 
Ni 3 Ga a = 140, b = 270, v F = 0.19, and £ = 0.556 Ms- 
Using the resulting values of £ each compound we ob- 
tain a magnetic moment of M = 0.3 /Us/cell for Ni 3 Al 
and a paramagnetic state with the renormalized suscep- 
tibility x(0,0) = l/a 2 = 6.8 x 10~ 5 emu/g for Ni 3 Ga, 
thus correcting the incorrect ordering of the magnetic 
tendencies of these two compounds and reproducing ex- 
tremely well the experimental numbers of M = 0.23 /is, 
x(0, 0) = 6.7 x 10~ 5 emu/g, respectively. This qualitative 
behavior is due to the different coefficient a, i.e., different 
q dependencies of Xo(<7, 0) at small q, which relates to the 
phase space available for soft fluctuations. 

Now we turn to the pressure dependence. The above 
results imply that beyond LDA fluctuations are already 
larger than the moments themselves at P = 0. In this 
regime, we may assume that the size of the beyond LDA 
fluctuations is only weakly pressure dependent. Then 
we can apply Eq. 2 to the data shown in Fig. 3 using 
£ = 0.47 as needed to match the P = value of M. This 
yields a value P c =10 GPa in quite good agreement with 
the experimental value, P c =8.1 GPa. [7] 

In conclusion, we address the LDA failure to describe 
the physics of magnetism in Ni 3 Al and Ni 3 Ga even qual- 
itatively. We identify the problem as neglect of spin fluc- 
tuations associated with the ferromagnetic quantum crit- 
ical point. These are stronger in Ni 3 Ga despite the fact 
that the latter has a larger density of states and is there- 
fore more magnetic in mean-field theories. The reason 
for the difference in the spin fluctuation spectra is in the 
q dependence of the non-interacting spin susceptibility. 
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FIG. 3. FSM calculations under hydrostatic pressures. 
Magnetic energy, defined as the energy relative to the 
non-spin-polarized result at the same volume, as a function 
of the moment and linear compression. Left and right panels 
correspond to Ni 3 Al and Ni 3 Ga, respectively. 
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TABLE I. Magnetic energy (see text), magnetic moment 
in /is/cell and N(Ef) in eV _1 on a per spin per formula unit 
basis. 





\AE\ (meV) 


M (calc.) 


M(expt.) 


N(E F ) 


Ni 3 Al 


10.3 


0.71 


0.23 


3.2 


Ni 3 Ga 


14.3 


0.79 


0.00 


3.4 
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